Fault-tolerant quantum computation with long-range correlated noise 
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We prove a new version of the quantum accuracy threshold theorem that applies to non- 
Markovian noise with algebraically decaying spatial correlations. We consider noise in a quantum 
computer arising from a perturbation that acts collectively on pairs of qubits and on the environment, 
and we show that an arbitrarily long quantum computation can be executed with high reliability in 
D spatial dimensions, if the perturbation is sufficiently weak and decays with the distance r between 
the qubits faster than 
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The theory of quantum fault tolerance [1] establishes 
that quantum information, if cleverly encoded, can be 
protected from damage and processed reliably with im- 
perfect equipment. One noteworthy achievement of this 
theory is the quantum threshold theorem [2-6], which as- 
serts that an arbitrarily long quantum computation can 
be executed with high reliability, provided that the noise 
afflicting the computer's hardware is weaker than a cer- 
tain critical value, the accuracy threshold. 

Early proofs of the threshold theorem [2-4] assumed 
that the noise is Markovian, that each quantum gate in 
the noisy circuit is a trace-preserving completely positive 
map that approximates the ideal gate. Each gate in this 
noise model can be realized as a unitary transformation 
that acts jointly on a set of the qubits in the computer 
(which we will call system qubits) and on the environ- 
ment (the bath variables), but where it is assumed that 
the bath has no memory — the state of the bath is re- 
freshed after every gate. In [7] the theorem was extended 
to a non-Markovian noise model of a restricted type. It 
was assumed that each qubit Qi has a local bath Ei , and 
that Ei and Ej interact only when a quantum gate acts 
on the corresponding pair of qubits. In [5], the noise 
model was relaxed further by allowing nonlocal interac- 
tions among the bath variables, but it was still assumed 
that a pair of system qubits interact directly only when 
a quantum gate acts on that pair. In this paper, we will 
show that the threshold theorem applies even if there are 
long-range interactions among the system qubits that are 
"always on." 

The noise model we will consider can be formulated in 
terms of a time-dependent Hamiltonian H that governs 
the joint evolution of the system and the bath. We may 
express H as 

H = Hs + Hb + Hsb , (1) 

where Hs is the time-dependent Hamiltonian of the sys- 
tem that realizes the ideal quantum circuit, Hb is an ar- 
bitrary Hamiltonian of the bath, and Hsb couples the 
system to the bath. In [5] it was assumed that the 
system-bath coupling could be expressed as 



Hsb — ^ HsB,a ; (2) 

a 

here a labels a circuit "location" during a particular time 
step of the computation — either a single qubit at a lo- 
cation where a one-qubit gate acts in the ideal circuit 
(which could be the identity gate in the case of a "rest- 
ing" qubit), or a pair of qubits at a location where a 
two-qubit gate acts on the pair in the ideal circuit. This 
coupling of the system qubits to one another is "short 
range" in the sense that two system qubits interact di- 
rectly with one another only if the ideal Hamiltonian Hs 
also couples those two data qubits during the same time 
step — hence we will refer to it as the short-range noise 
model. For short-range noise it was shown in [5] that 
there is a positive constant Eq such that an arbitrarily 
long ideal quantum computation can be accurately sim- 
ulated using the noisy gates (with reasonable overhead) 
provided that at each time and at each location 

II HsB,a II to <eo , (3) 

where to is the time needed to execute a quantum gate 
and [I'll denotes the sup operator norm. 

In this paper we will consider a system-bath coupling 
of the form 

Hsb = E • (4) 

where the sum is over all pairs of system qubits; Hij acts 
collectively on the pair of qubits < ij > and on the bath. 
(By convention, Hij = Hji — H.^ijy.) Since any pair 
of system qubits can interact, we call this the long-range 
noise model. We will show that for long-range noise there 
is a positive constant rjQ such that an arbitrarily long ideal 
quantum computation can be accurately and efficiently 
simulated provided that for each qubit i and at each time 

J2 II H,, II to < % . (5) 

j 

We note that if the interaction between system qubits 
decays algebraically so that 
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\\Hii\\<5/\i-jW (6) 

where |i — j| denotes the distance between qubits i and j, 
then the sum in eq. (5) converges in the hmit of an infinite 
f-dimensional system provided that z > D. Thus, our 
result shows that fault-tolerant quantum computation is 
possible for 6 sufficiently small. 

As in [7,5], we will analyze the noisy circuit using a 
"fine-grained fault-path expansion." Wc divide the work- 
ing period to of a gate into to/A:^ 1 intervals each of 
width A; for short-range noise the time evolution opera- 
tor for the time interval {t, t + A) can then be expressed 
as 

U{t + A, i)« e-"^" 

« e-^^Hs^-^AHs -Q ^j^^ _ ,AHsB,a) (7) 
a 

since for A sufficiently small terms higher order in A 

can be safely neglected. The action C/circuit on the sys- 
tem and bath of a noisy quantum circuit with altogether 
N time steps can be expressed as a product of Nto/A 
such cvohition operators. The fine-grained fault-path ex- 
pansion of the noisy circuit is constructed by expressing 
^^circuit as a sum of monomials, where in each term cither 
IsB or —iAH SB, a is inserted at each "micro- location" 
with duration A. A micro- location where —iAHsB,a 
is inserted is called a "micro-fault." A noisy gate or 
"macro-location" in the circuit consists of to/ A consec- 
utive micro-locations, and the macro-location simulates 
the ideal gate perfectly if it contains no micro- faults. 
Therefore, we say that a macro-location is faulty (is a 
"macro-fault") if it contains one or more micro-faults. 
By grouping together terms in the fine-grained fault-path 
expansion we obtain a coarse-grained fault-path expan- 
sion such that each term in the coarse-grained expan- 
sion is the sum of all terms in the fine-grained expansion 
with macro-faults occuring at a specified set of macro- 
locations. 

The strength of the noise in the circuit can be quanti- 
fied by the norm of the sum E{Ti) over all the terms in 
the fine-grained expansion such that a particular spec- 
ified macro-location li is faulty (where no restriction 
is placed on the noise at other macro- locations). We 
can organize this sum by noting that there must be an 
earliest micro-fault in the macro-location which can oc- 
cur at any of the to /A micro- locations; summing over 
all fine-grained fault paths with a fixed earliest micro- 
fault is equivalent to inserting Is b at the micro- locations 
prior to the earliest micro-fault and inserting e-^^^sB.a 
at the later micro-locations. Since || e-^^^sn.a ||_ i 
and the norm of a product is bounded above by a prod- 
uct of norms, the norm of the sum over all fine-grained 
paths with a fixed earliest micro-fault is no larger than 
A II HsB,a II- Now summing over the to /A possible 
micro- locations for the earliest micro- fault, and noting 
that the norm of a sum is bounded above by a sum of 
norms, we have 



II E{Ii) II < £ = Max II HsB,a II to , (8) 

where the maximum is with respect to all a and all times. 
By similar reasoning, if E{Ir) denotes the sum over all 
fine-grained fault paths such that faults occur at a par- 
ticular set Xr of r macro-locations (where no restriction 
is placed on the noise at other macro-locations), then 

II E{Xr)) II < 8^ . (9) 

If cq. (9) is satisfied, we say that the noise is "local 
with strength £." Note that e is determined only by 
the system-bath coupling Hsb', no restriction need be 
imposed on the bath Hamiltonian Hb- 

Thus we see that the short-range noise model obeys the 
local noise condition eq. (9). For a quantum computer 
subjected to local noise, the effectiveness of a recrusive 
fault-tolerant simulation of an ideal circuit was analyzed 
in [5]. In a level- 1 simulation, each gate of an ideal cir- 
cuit is replaced by a 1-gadget, a circuit that simulates 
the gate acting on quantum information protected by a 
quantum error-correcting code. In a level-fc simulation, 
each gate of an ideal circuit is replaced by a fc-gadget, 
which is constructed by replacing each gate in a (fc— 1)- 
gadgct by a 1-gadget. If the code can correct t errors, 
then a 1-gadget is said to be bad if it contains t + 1 or 
more faulty macro-locations, and a fc-gadget is said to be 
bad if it contains t + 1 or more bad (fc— l)-locations. 

If the 1-gadgets are properly constructed, then it can 
be shown that for a level-1 simulation, the norm of the 
sum E{2r-^^) over fault-paths such that the 1-gadgets are 
bad at a specified set Jr^^ of r level-1 locations satisfies 

II E{4'^) II < (^(i))' , (10) 

where 

here A is the maximal number of macro-locations con- 
tained in any 1-gadget, and C is a constant of order 1. 
Furthermore, by an inductive argument we find that for 
a level-A; simulation, the norm of the sum over 
fault-paths such that the fc-gadgets are bad at a specified 
set Jr'^^ of r level-A: locations satisfies 

II ^(XW) II < (eC^))' , (12) 

where 

e^'^=c[^^,){e^'-'^f'=eoie/sor'^' ' ^''^ 
and 
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Thus, for e < Eq, the norm drops steeply as the level 
k increases, and it can be shown that an ideal circuit 
with L gates can be simulated with fixed accuracy and 
an overhead cost that is poly logarithmic in L. This is 
the quantum threshold theorem for local noise [5]. 

To establish a threshold theorem for the long-range 
noise model, then, it will suffice to show that the local 
noise condition cq. (9) holds for the long-range noise, 
just as it does for short-range noise. For now, suppose 
that all of the gates performed at the r specified faulty 
macro-locations are single-qubit gates. If these macro- 
locations all occur at different times, then our earlier 
analysis of short-range noise needs little modification. 
Now we express the time evolution operator for time in- 
terval {t, t + A) as 

U{t + A, i) « e-'^"'e-'^"'' II ^^Hij) (15) 

<ii> 

and again expand ?7circuit to generate a fine-grained 
fault-path expansion. Considering a fixed faulty macro- 
location Xi where the ideal gate acts on qubit i, again 
there vaxist be an earliest micro-faiilt that acts on qubit 
i in that macro-location, which can be at any of Iq/ A 

intervals; we insert — iA Hij^ at that earliest micro- 
fault, I at each earlier micro-location, and q-'^'^^sb q\i 
later time intervals. Then summing over the position of 
the earliest micro-fault we find 

II £(Ji) II < ?? = Max II Hij II j to , (16) 

where the maximum is with respect to all i and all times. 
Similarly, if macro-faults occur at a set of Ir of r macro- 
locations, where all the macro-locations are at different 
time steps, then by summing independently over the ear- 
liest micro-fault at each macro-location, we find 

II E{Ir) II < V' ■ (17) 

But what if many macro-faults occur during the same 
working period? The long-range noise model allows faults 
occuring in the same time step to be correlated. The 
worst case is when the r specified faulty macro-locations 
in the set 2^ are all in the same time step. In that case, if 
i and j are two of the qubits acted upon by the r faults, 
then Hij can account for the two of the faults. 

We again organize the sum over fine-grained fault 
paths by considering the earliest micro-fault that strikes 
each of the r qubits. For a particular fine-grained fault 
path, if i, j e and Hij is the first perturbation to act 
on both of qubits i and j, then we say that the pair ij 
is "contracted." If Hij is the first perturbation to act 
on qubit i, where either j Ir or j G Xj. has already 
been struck earlier by another perturbation, then we say 
that qubit i is "uncontracted." By definition, then, the 
contracted pairs are disjoint. 



We group the fine-grained fault paths into families, 
where in each family the time of the earliest micro-fault 
acting on each qubit in X,. is specified, as well as the pair 
of qubits acted upon by that fault. If is clear that the 
norm of the sum of fine-grained fault paths in a fam- 
ily can be bounded by a product — with a factor of 
A II Hij II arising from a fault that acts on the pair 
< i] >; the only problem is to count the families. 

In each family, a specified set of fc pairs are contracted, 
where k ranges from to r/2 for r even, or from to 
(r — l)/2 for r odd. Each contraction can occur in any 
of the io/A time steps, which converts the factor of A to 
tQ. For an uncontracted qubit i, we sum over the qubit 
2 it could be paired with by the perturbation Hij — the 
sum runs over both j G and j ^X^. For j ^ X^j this 
micro-fault can occur at any time, but for j € X^ that 
might not be true. If we slide Hij forward in time to 
before the first micro-fault acting on qubit i, then we are 
considering a family in which i and j are contracted. But 
that is all right; if we allow Hij to occur at any time, we 
are only overestimating the norm of the sum over fault 
paths; so here again we have an upper bound by convert- 
ing A to to J and summing over j gives a factor of r\ for 
each uncontracted qubit. 

Let's suppose that r = 2n is even. The norm Sk of the 
sum over fine-grained fault paths with k contractions can 
be bounded above: 

Sk < Yi n (II II *o) 

contractions i=l 

^FH <""'''* 

(Note that the sum of || Hij \\ over all i and j, raised to 
the power k, contains a term bounding the contribution 
from each contraction 2^^^ times, and also contains other 
positive terms.) Thus we have 

(W^ Vk) ■ 

To bound the sum over k we use 




and n 4- 1 < e" to obtain 

n 

II E{Xr) 11= ^ 5fe < (n + l)e"e"/^?7" 

fe=0 

< (e^+Ve^)" = ^ (21) 

where 
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V = 6^+1/2'= 0? • (22) 

Thus, wc have shown that, even when aU of the r — 2n 
specified faults occur during the same working period, 
the norm of the sum over fault paths is bounded by tj"', 
where r}' = 0{^). This bound also applies when r is 
odd. 

The same bound applies to two-qubit gates, except 
that we need to replace r] by 2r]. A two qubit gate has a 
fault if the perturbation hits either one of the two qubits 
(and after one of the qubits is damaged, we don't care 
what happens to the other one). The first perturbation 
to damage the pair of qubits i and j might act on the 
pair im (m ^ j), on the pair mj (m =/= i), or on the pair 
ij. So the norm of the sum over fault paths that damage 
one particular two-qubit macro-location can be bounded 
above by 

I E II II + E II II + II \\]to<2v. (23) 

When we consider summing over the fault paths with 
many faults in the same time step, where some of the 
faults occur at two-qubit macro-locations, now the "un- 
contracted" contributions include the case where a sin- 
gle perturbation damages both of the qubits at a single 
macro-location; if on the other hand a two-qubit macro- 
location is contracted with another macro-loeation, the 
perturbation can act on either qubit at the two-qubit 
macro-location. Both types of contributions are taken 
into account by replacing 77 by 2ri. 

The conclusion is that "long-range noise" is really a 
type of "local noise" obeying eq. (9), where 

e = ei+i/2-y2^ . (24) 

therefore the quantum accuracy threshold theorem 
proven in [5] is applicable. Note that we can use simi- 
lar reasoning for many-body interactions. For example, if 
the system-bath interaction is a sum of three-body terms, 
then 

r?= ^MaxE II i^ijfc llj io , (25) 

and e = 0{-q^/^). 

Actually, even the "short-range" noise model includes 
long-range spatial correlations among the faults. If 
each system qubit couples individually to the bath with 
strength A. but the self-interactions of the bath may be 
nonlocal, then a coupling between a pair of system qubits 
of order can be induced in second-order perturbation 
theory, and p-body terms of order A'' are generated in 
higher order. Nevertheless, the "long-range" noise model 
is more general, because there are interactions among the 
system qubits that are allowed in the long-range model 
but cannot be simulated by a local coupling of the sys- 
tem to the bath. Furthermore, by working directly with a 



system-bath Hamiltonian that includes long-range inter- 
actions among the system qubits, we can see more easily 
how the rate of decay of the interactions with distance 
enters into the analysis, as in eq. (6). 

Using methods described in [5], we can obtain a lower 
bound on the threshold noise strength £0 ~ 10~^, or 
?7o » 10~^° for the case of two-body interactions. This 
estimate is discouraging, but it is probably much too pes- 
simistic, for several reasons. In the analysis we allow the 
bad fault paths to add together with a common phase, 
so that they interfere constructively. Under a more plau- 
sible scenario in which the phases of distinct fault paths 
are only weakly correlated, the norm of the sum is sig- 
nificantly smaller, so that much stronger noise can be 
tolerated in practice. Furthermore, by following sound 
engineering principles, for example by spatially separat- 
ing qubits belonging to the same code block, we can im- 
prove the effectiveness of fault-tolerant protocols. We 
also note that more clever gadget design could signifi- 
cantly improve estimates of the threshold as suggested 
for example in [8] , and in any case the actual value of the 
threshold is likely to be notably higher than can be es- 
tablished by rigorous arguments. The noteworthy point 
is that we have shown that sufficiently weak long-range 
interactions among the qubits (and interactions with a 
nonlocal bath that can store information for an indefi- 
nitely long time) are compatible in principle with fault- 
tolerant quantum computation. 
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